A new theoretical approach to Interacting Dimers 
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Recent numerical results on classical dimers with aligning interactions have been theoretically 
explained via a Coulomb Gas representation of the height random variable. Here we propose a 
different, more direct, theoretical viewpoint, the Interacting Fermions Picture, which provides a 
better account of the large-distance asymptotic formulas of certain correlations. 
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I. INTRODUCTION 

The lattice model of hard-core close-packed dimers is 
among the most fundamental systems in two-dimensional 
Statistical Mechanics. Not only it is exactly solvable in 
a broad sense (the free energy and many correlations can 
be computed, [i- r >]) but it also provides, through equiv- 
alences, the exact solutions of several other models, in- 
cluding the nearest neighbor Ising model, [1, 6], and some 
vertex models at the so called free fermion point, [7, 8] . 

Recently, especially in connection with a problem of 
Quantum Statistical Mechanics, [9], several authors have 
been studying the (classical) dimer model on a square 
lattice, modified with a nearest neighbor aligning inter- 
action. We will call it interacting dimer model, (IDM). 
Since a general exact solution for the IDM has not been 
found so far, our knowledge of the properties of this sys- 
tem - quite a deep knowledge, in fact - rests entirely on 
the numerical analysis of [10-13]; and on their theoretical 
interpretation via the Coulomb Gas approach, (CGA). 

The CGA, [14-16] (see also [17] 's review), applies ev- 
ery time in the model there is a natural height random 
variable; and it is based of the postulate that, in the scal- 
ing limit, height variable correlations are equal to charge 
correlations of a free boson field. From a practical view- 
point, the CGA has been very successful in describing 
the critical exponents of virtually every critical model of 
two-dimensional, lattice Statistical Mechanics. Though, 
the conjectured scaling limit of the height variable has 
been difficult to substantiate (the best result in this di- 
rection is still nowadays [18]). Besides, in the case of the 
IDM, the CGA has not provided so far a clear account 
of some fine details, such as the staggered prefactors of 
the correlations (see below). 

This article proposes a different theoretical approach 
to the IDM, that can be called Interacting Fermions Pic- 
ture (IFP). This approach is not new in Physics. It 
has been a standard tool in condensed matter theory for 
studying the 1+1-dimcnsional quantum models (see re- 
views [19, 20]). In classical 2-dimensional Statistical Me- 
chanics, it was first employed in [21] to demonstrate the 
universality of the nearest neighbors Ising model under 
small, "solvability breaking" , perturbations; and recently 
it was used to study the weak-universal properties of the 
Eight- Vertex and the Ashkin- Teller models, [22, 23]. In 
regard to the IDM, fermion viewpoints have already been 
employed in [24, 25]; however, the full fledged IFP ap- 
plication to the IDM that is given here appears to be 



completely new. The method is made of two steps: i) 
the IDM is re-casted into a lattice fermion field with a 
self-interactions; ii) the scaling limit of the lattice field 
is proved to be the Thirring model, which is interact- 
ing as well, but which is also exactly solvable. The IFP 
solves some difficulties of the CGA because: a) the con- 
vergence, in the scaling limit, of lattice correlations to 
the Thirring model correlations naturally emerges from 
renormalization group ideas; b) it makes clear the origin 
of the staggering prefactors of the dimer correlations. Be- 
sides, via bosonization of the Thirring model, one could 
always recover the free boson description that is typical 
of the CGA, without ever introducing the height variable 
of the lattice model. 

In our view, the IFP is a most natural and transparent 
way to describe the IDM. 



II. DEFINITIONS AND RESULTS 

Consider a finite box A of the infinite square lattice. A 
dimer configuration. u>, is a collection of dimers covering 
the edges of A with the constraint that every vertex of 
A is covered by one, and only one, dimer. The partition 
function of the interacting dimers model (IDM) is 



*(A) = £' 



AS 



, v{d,d'\ 



(i) 



where: A is a real parameter; v(d, d') is a two body dimer 
interaction; the first sum is over all the dimer config- 
urations; the second sum is over any pair of dimers in 
the configuration cj. In [11] \v{d,d') is a nearest neigh- 
bor aligning interaction. In this work we assume more 
in general that v(d,d') is zero unless d and d' are both 
horizontal or both vertical; that it is invariant under 
7r/2-rotations and under lattice translations; and that 
\v(d, d')\ has exponential decay in the distance between 
d and d! . The non-interacting, exactly solvable, dimer 
model is the case A = 0, [1—5]. 

Our main result is the evaluation of correlation crit- 
ical exponents of local bulk observables, for small |A|. 
A natural observable to consider is the dimer occupancy 
Vdito), which is equal to 1 if the dimer d is present in u>, 
and it is zero otherwise. Consider the horizontal dimers 
d = {0, eo} and d' = {x, x + eg}, for eo = (1,0) and 



x = (xq,xi). The IFP provides the large- |x| formula 



Xq+X! 



{VdVd') - (vd){vd') ~ (-1) 



+ (-l) a °c_ 



(org + x\Y- 



(2) 



for critical exponent k_ = 1 + 0(Ai>) (and for Xv- 
dependent c and c_ ) . Notice that according to the sign of 
A cither the former or the latter term in (2) is the leading 
one at large |x|. (2) coincides with the exact solution for 
A = 0: see (7.12) and (7.20) of [3]; and it is in agreement 
with the numerical simulations for A > (and smaller 
that a threshold value): see (51) and (52) of [11]. 

The critical exponents n_ is non-universal, because it 
does depend upon \v. What is expected to be univer- 
sal, instead, is the relationship among critical exponents 
of different observables. It is instructive to study a sec- 
ond observable, then. The authors of [11] considered the 
monomer correlation. Wc could describe this correlation 
in the IFP as well; however, being the monomer corre- 
lation equivalent to a non-local fermion correlation, the 
derivation of the scaling limit in the IFP, at the present 
time, is not more transparent than in the CGA. 

We consider instead a different observable, the diago- 
nal dimer. It consists in a pair of monomers in positions 
{x, x + e}, where e = (1,1); see Fig.l. Since such a 
dimer is not allowed in the hard-core, close-packed con- 
figurations uj, we define its "correlation" as it is done for 
the monomer observable, i.e. in terms of lattice defects: 



Z X (A - d) 
( ^ )= A lim o Z X {A) {VdVd ' 



lim 



Z x (A-(dUd')) 
Zx(A) 



For d = {0,e} and d' = {x, x + e}, the IFP gives the 
large- 1 x| formula 



(-l) x ° - (-1) X1 



(4 



(3) 



for a new critical exponent n + = 1 + O(Xv) (and for a 
Aw-dependent c+). The universal formula that relates /c+ 
to k_ is peculiar of the models with central charge c = 1 
and was originally discovered by Kadanoff, [26]: 



= 1 



(4) 



In the next section we introduce the fermion representa- 
tion of the IDM and derive (2), (3) and (4). In fact, a 
power series expansion in A for k_ and k+ could also be 
provided, but we will not do it here (see [23]). We re- 
mark that the CGA was successfully used to justify the 
appearance of critical exponents k_ and k+ which satisfy 
(4) (see pt.l and pt.5 on page 7 of [27]); however, in that 
approach, the derivations of the staggering prefactors of 
(2) and (3), as well as of the power series formulas for k_ 
and k+, are not available. 



FIG. 1: Example of a diagonal dimer. Its presence in the 
graph generate a new face with: 6 vertical arrows, 3 of which 
are clockwise oriented; 6 horizontal arrows, 3 of which are 
clockwise oriented. The total oddness of this face is then 
i 3 ( — i) 3 l 3 ( — l) 3 = —1 and so, as opposed to the case of the 
monomer defect, [3], the fermion representation of the diago- 
nal dimer remains local. 



III. INTERACTING FERMIONS PICTURE 

When A = the dimer model is equivalent to a lattice 
fermion field without interaction. Namely 



Z (A) = JDipexp{-^Yl K *,y^y} 



(5) 



where: {ip x ■ x G A} are Grassmann variables and Dip 
indicates the integration w.r.t. all of them; K x y is the 
Kasteleyn matrix that can be chosen to be such that 



/ ,-^x,y 



x+creo 



+ ^x+crei ) 



er=±l 



with eo = (1, 0) and ei = (0, 1). (5) is the partition func- 
tion of a free Majorana fermion field, i.e. a Grassmann- 
valued Gaussian field with moment generator 

where: {n x : x e Z 2 } are other Grassmann variables; S, 
the covariance, is the inverse Kasteleyn matrix 



^ipoxa+ipixt 

2n J_ n 2tt isinpo — sinpi 

The Fourier transform of S is singular at four Fermi mo- 
menta: p+,o = (0,0), p+,i = (n, 7r), P-,o = (0, 7r) and 
p_,i = (tt, 0). Therefore, in view of the scaling limit, it 
is convenient to decompose 



E 

u> = ± 
s = 0,l 



where ip x . UtS are four independent Majorana fields with 
large- |x| covariances 

1 [dpo [dpi jPovo+ipixi 

{lpx,u,s1p0,w,s) 



2 J 2tt 
1 1 

Air xq + iojX\ 



2tt 



ip - bjpt 



(6) 
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This decomposition already appeared in [24]. Anyways 
here we are preparing for the application to the interact- 
ing case, therefore we also introduce the Dirac spinors 
i't = W£+,^x,-) and V^x = (^x,+ ,^x ,-) T for 



,+ _ V'x.c^O + #x,id,l , _ 1px,u,0 - ii>x,u, 
Vx,lu — rx.u — 



y/2 rX '" V2 

with translational invariant covariances 



i (7) 



i 



(8) 

47T Xq + ILUXl 

If we now let A ^ 0, by power series expansion in A one 
can verify (by cluster expansion) that (1) becomes 

Z X (A) = flhl> CX P { - 1 Y, K x,yMy + V(XV, V-)} 
J x,y 

(9) 

where V(\v,?p) is an even polynomial in the tp. In the 
next section, by a Renormalization Group argument, we 
will explain why, in the evaluation of the large distance 
decay of the correlations, it is correct to replace the in- 
teracting fcrmion field (9) with the Thirring model. As- 
suming for the moment this crucial fact, we now compute 
the correlation critical exponents. It is not difficult to see 
that the dimer correlation in the l.h.s. of (2) becomes, 
in terms of the Dirac fermions (7) (and up to terms with 
faster decays) 

(V'oV'eoV'xV'x+eo) - (V'oV'eoXV'xV'x+eo) 

~ 4(-i) a °- HBl £ «^o,.;Ox,.) T 



+ 4(-l)*° Y «>o,-c; i>t- u ^) z 



(10) 



where the label T indicates a truncated correlations. In 
the same way, the diagonal dimer correlation in the l.h.s. 
of (3) becomes 

(^oV'eV'xV'x+e) - (ipQlpe) (ipxi'x+e) 

~ 8 - (-1H « + ^o + _;^x,-^x, + ) T (ii) 

The last step is to look up the exact solution for the 
massless Thirring model correlations, [28-30] (see also 
the review [31]). We find 

(^(0)^(0); Vi(*)^(x)) T 



2\2 



(^(0)Vl(0); ^_(x)^+(x)) T = 
(^(0)^_ u (0);^ u (x)^(x))' 



{xl+x\) 
c_ 

(xl + x\Y+ 
c+ 

(xi + x \y- 



(12) 



for critical exponents 



i + 4 21 

i At 

1 4?r 



i — An 

1 + At 



where Ay = O(Aw) is a parameter of the Thirring model. 
The derivation of (2), (3), (4) is complete. 



IV. RG ANALYSIS 

We follow the Wilson's RG scheme in the version due 
to Gallavotti, [32]. Integrating out the large momentum 
scales, one obtains an effective interaction 



E 



E 



dki ■ ■ ■ dk 2n ^k u uj uSl ■ ■ ■ ^k 2 „, 



' <H E k J + E PtJ i' s 0^ 2l ^ k2 +PtJ 2,S2;- • • > k 2« + Pw 2 „,s 2n ) 

It is important to take into account some symmetries. 
For R(ko,ki) = (ki,—ko), i?(fco,fci) = (fci,fco) and 
r(fco, fci) = (fco, fei + 7r), we find 

W2m(rk 2 , . . . , rk 2m ) = (-z) m w 2m (^k 2 , . . . , $k 2m ) 

w 2m (i?k 2 , . . . , i?k 2m ) = (-z) m w 2m (k 2 , . . . ,k 2m ) 

w 2m (k 2 , . . . ,k 2m )* = i m to 2m (i?k 2 , • • ■ ,#k 2m ) . (13) 

From power counting, there are two possible local, 
marginal terms: a quartic term, that requires the renor- 
malization of the coupling constant A; a quadratic term, 
responsible for a field renormalization. By (13) they are: 



4u5 4 (p +! i, p_ i0 , P-,i) E ^i+V'x.+V'x"-^ 



and 
2 



. dw 2 , \ dw 2 . 



r: 

X,CJ 



(14) 



(15) 



where is the Fourier transform of iko — u)k\. Again 
by (13), the prefactors in (14) and (15) are real. In- 
stead, there are no local, relevant terms: the only pos- 
sible one, a quadratic terms without derivatives, i.e. a 
mass term, cannot be generated by conservation of the 
total momentum. These facts imply that (9) (with pa- 
rameter A) equals the massless Thirring model with pa- 
rameter Ay, up to terms which are irrelevant and thus 
cannot break down (4) (although they do determine the 
functional relationship between A and Ay). 



V. CONCLUSION 

We have showed that the IFP provides a detailed the- 
oretical account of the numerical findings of [11] -plus 
some new predictions- in the case of weak interaction. 
The method should also work, with possibly different 
outcomes, for triangular and hexagonal lattices; and for 
two or more interacting copies of dimer models. Besides, 
the IFP should be applicable to the Six- Vertex model, 
which is equivalent to dimers on a square lattice with a 
staggered interaction. Including the results on Ashkin- 
Tellcr, Eight- Vertex and XYZ quantum chain, [23], the 
IFP seems quite an effective way for dealing with two 
dimensional, lattice, critical models with central charge 
c = 1. It might be possible that the IFP be applicable to 
c < 1 models (for example through ideas in [15]). 
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VI. ADDED NOTE: STRONG INTERACTION 

In the opposite case of strong dimer interaction, i.e. 
|A| ^> 1, the numerical findings of [If] indicate the exis- 
tence of four different "columnar" phases. A theoretical 
justification of this case is also possible; however, as op- 
posed to the weak interaction case, the outcome crucially 
depends on the choice of the interaction v(d, df): for def- 
initeness, we only discuss here the choice in [11], which 
assign an energy — A < per each plaquette with one of 
the two dimer arrangements in Fig. 2. 

Decorate a dimer configuration with wiggly lines as 
indicated in the same figure. Declare a plaquette to be 
"good" (or "ground state plaquette" ) if it contains 2 wig- 
gly lines; otherwise, declare the plaquette to be "bad" (or 
"excitation plaquettes" ) . Note that: a) dimer configura- 
tions on nearest neighbor good plaquettes must corre- 
spond to the same columnar ground state; b) the prob- 
ability of a bad plaquettes occurrence is dumped by a 
factor e~~, or smaller, per plaquette. Therefore, for A 
positive and large, one can apply the standard Peierls' 



argument to prove the existence of four different phases 
(see [33, 34] or the review [35]). 





FIG. 2: For each plaquette containing two facing dimers 
draw four wiggly lines. An equivalent, but completely local 
way of evaluating the total energy of a dimer configuration is 
by assigning an energy of — h per wiggly line. 
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